Designing highly efficient interlocking interactions in anisotropic active particles

Cluster formation of microscopic swimmers is key to the formation of biofilms and colonies, efficient motion and nutrient uptake, but, in the absence of other interactions, requires high swimmer concentrations to occur. Here we experimentally and numerically show that cluster formation can be dramatically enhanced by an anisotropic swimmer shape. We analyze a class of model microswimmers with a shape that can be continuously tuned from spherical to bent and straight rods. In all cases, clustering can be described by Michaelis-Menten kinetics governed by a single scaling parameter that depends on particle density and shape only. We rationalize these shape-dependent dynamics from the interplay between interlocking probability and cluster stability. The bent rod shape promotes assembly in an interlocking fashion even at vanishingly low particle densities and we identify the most efficient shape to be a semicircle. Our work provides key insights into how shape can be used to rationally design out-of-equilibrium self-organization, key to creating active functional materials and processes that require two-component assembly with high fidelity.

The motility inherent to synthetic [1][2][3] and living [4][5][6] microswimmers is crucial for their collective behavior.Besides pushing these systems out of equilibrium, it is at the heart of the occurrence of trains 7,8 , flocks 3,9 , vortices 9,10 , and clusters 1,2,[4][5][6]11 . Cluter formation occurs in a range of active agents, from self-propelled colloids forming living active crystals to rotating or moving clusters of populations of starfish embryos 6 , or motile bacteria like Thiovulum majus 5 and Myxococcus xanthus mutants 4 .There, it is thought to be important for the formation of colonies and biofilms, efficient motion, and nutrient uptake 12 .
Simulations of self-propelled agents have identified that motility is a minimal requirement for cluster formation, and that it occurs even in systems of repulsive agents at sufficiently high densities.This socalled motility induced phase separation (MIPS) 11,[13][14][15][16] originates from a slowing-down of the active agents upon collision.Systems of selfpropelled disks at high particle densities phase separate above a critical density of 40% 17 , forming a single, globally disordered, macroscopic MIPS-aggregate that experiences diffusive motion.In experiments with synthetic microswimmers, however, cluster formation usually sets in already at several percent due to the presence of additional attractive interactions, and multiple disordered and dynamic aggregates are found across the sample 1,2,18 .
In contrast to synthetic active particles, biological microswimmers often have non-spherical shapes [19][20][21][22] and simulations have predicted that their collective behavior is strongly influenced by shape.For example, longitudinally-propelled rods self-organize into polar bands instead of clusters as a consequence of the lateral association promoted by their elongated shape [23][24][25][26] .Clusters of anisotropic active particles may not only display Brownian diffusion, but also directed and spinning motion 27,28 .In addition, the critical density at which cluster form has been proposed to depend on shape because it strongly influences inter-particle slowing down 28 .For active hexagons, for example, efficient deceleration upon collision results in nucleation of many small clusters at far lower particle densities than would be expected for spheres 28 .
Despite these predictions, there is little experimental work due to the limited availability of active particles with an anisotropic shape.Experiments with transversely-propelled rods revealed that stable doublets already form at particle densities of ≈ 1% 29 and clustering is favored with increasing aspect ratio of the rods 29,30 .Tori, both horizontally and vertically oriented, were observed to form dynamic unstable clusters 31 and standing disks were found to cluster at about 10% 32 .While more complex shapes have been prepared 33,34 it is often challenging to selectively generate sufficient quantities to test their phase behavior and typically not possible to gradually tune their shape to find optimal clustering conditions.
Here, we use 3D micro-fabrication to create active particles with a shape that is continuously tunable between a sphere and a rod, i.e. bent rods, and study their collective behavior at low particle densities.We find that clustering occurs at extremely low surface area fractions (below 0.022%) due to a highly efficient nucleation process that relies on interlocking between the particles.Complementing our experiments with simulations and an analytical model, we find that the self-organization at low particle densities can be captured by a Michaelis-Menten kinetics and thus can be characterized by a single scaling parameter that depends on the particle density and shape only.We demonstrate that the efficiency of the self-organization process is strongly influenced by the concave shape, which affects the interplay between interlocking and stability, and show that the combination of directed motion and steric interactions leads to motility induced clustering.Our insights provide a generic understanding of how phase separation occurs in a whole class of anisotropic particles from spheres to rods as well as a strategy to precisely control the stability of active particles through the shape of their interaction site, thus enabling the rational design of their assembly pathways into functional active materials.

Results
We exploit 3D microprinting based on two-photon polymerization to create active particles with a tunable anisotropic shape, see the Methods section for more details on particle fabrication.We print bent rods with varying opening angle α and constant cross-section L = 10 μm as this shape interpolates smoothly between a sphere and a rod, see Fig. 1a, b 35 .We render the particles active by sputter coating them on their convex side with a 5 nm thick Pt/Pd (80/20) layer and dispersing them in aqueous hydrogen peroxide (H 2 O 2 ) solution.Their self-propulsion is driven by solute gradients generated through a catalytic decomposition of H 2 O 2 at the Pt/Pd cap 36 .Depending on the fuel concentration, the bent rods either swim with their concave (1% H 2 O 2 ) or convex side (5% H 2 O 2 ) leading, see Supplementary Videos 1-3.Due to their size, the motion of these active crescents shows long persistence lengths, see Fig. 1c.

Collective behavior of active bent rods
We start by examining the collective behavior of active bent rods that swim with their concave-side leading, have an opening angle of α = 180°and an average velocity of 〈v〉 180 = 0.78 ± 0.08 μm s −1 .Small rotating clusters and pairs of particles quickly form after sedimentation, making our system effectively 2D, see Fig. 2a, b and Supplementary Videos 4-6.Simulations using a minimal model that only contains two ingredients, i.e., self-propelling bent rods and steric interactions between particles 37 , also show the quick formation of stable pairs and rotating clusters (Fig. 2c and Supporting Video 13).Clearly, the concave shape of the active particles is crucial as it promotes interlocking and stability against breakup.
Remarkably, stable interlocking occurs already for only two particles upon a head-on collision.The translational motion of the individual particles is then transformed into a stationary rotation of the pair, see Fig. 2d.The direction of rotation is determined by the torque created by the initial off-center alignment of the particles and remains stable over the duration of our experiment, with a typical angular velocity of ≈ 0.13 ± 0.01 rad s −1 .This angular velocity is in line with the active forces at play, which are in the range of F active = 0.04-0.1 pN (see SI).Not all collisions, however, will lead to the formation of clusters.The success depends on the probability for interlocking upon collision.For bend rods with an opening angle of α = 180°most collisions are successful: Besides direct interlocking upon collision, interlocking may occur after a short delay, i.e. when particles initially meet with their terminal parts and form a metastable ring-like structure.Any subsequent sliding induced by noise will result in interlocking, only with a delay.The quick formation and long stability of crescent-pairs seen in both experiments and simulations is very different from the behavior of pairs of active spheres, which tend to slide pass each other as soon as their velocities depart from a perfectly anti-aligned configuration.While cluster formation of spheres requires at least three, but typically more, particles 1,2,18 , the crescent shape favors configurations in which the end of one crescent is locked at the center of another particle and thus already stabilizes clusters of two particles.
The stable, rotating pairs subsequently act as nucleation points for larger clusters.Notably, clusters only grow through addition of single crescents and not by cluster merging, because of a lack of translational motion once the crescents form a pair.The long lifetime of the pairs enhances nucleation, similar to active polygons studied in simulations 28 .But unlike active polygons, active crescents can interlock as their shape features a cavity.Interlocking has so far only been reported for passive colloidal particles 38,39 or active systems sometimes in combinations with passive agents, neither of which were, however, consisting of only self-propelled and anisotropic colloidal agents [40][41][42] .Here, the interlocking makes the pairwise self-assembly stable against break-up, even upon a collision with a free particle, which is often absorbed into the cluster owing to its concave shape.We quantify the evolution of the cluster distribution over time and show the result for concave-side leading active bent rods (at α = 180°, 1% H 2 O 2 ) in Fig. 2e.Already after 10 min about 17% of crescents have formed pairs, despite our experiments taking place at a very low particle density of only ϕ s = 0.261%, which corresponds to only about 104 particles/mm 2 .Over the next 80 min the amount of single particles drops to 16%, more than half the particles are paired, while three-and four-particle clusters form in smaller amounts (20% and 9%, respectively).
This efficient nucleation and growth process leads to a quick assembly of the active bent rods into clusters.We find that the time evolution of the fraction of particles that are part of a cluster, N c , relative to the total number of particles in the system, N, quickly grows such that the majority of particles are part of a cluster after only about 30 min, i.e.N c /N > 0.5 (Fig. 2f), and almost all particles are part of a cluster after 90 min (Fig. 2e) despite the very low particle density of ϕ s = 0.261%.
To better understand the observed clustering process we consider a balls-into-bins model, where particle clustering is modelled as the presence of multiple particles (balls) inside a finite set of containers (bins) after random assignment (see SI).For each time step, all single particles are randomly assigned a new bin.A n-particle cluster corresponds to a bin containing particles.Simulating many time steps, we find a N c /N 0 curve as shown in Fig. 2g, which is qualitatively similar to those of the experiments.
To obtain an analytical expression, we consider only twoparticle-cluster formation for simplicity, because dimers are e Time evolution of the probability P(n) to find a crescent in a cluster of size n (ϕ s = 0.261%, measured after 90 min).The average particle velocity does not change significantly even after long times (> 5h, see SI). f Time evolution of fraction of particles that are clustered, N c /N, for the data shown in (e).Eq. ( 1) was used as a fit with τ N = 25 ± 4 min.g Balls-into-bins model with N 0 = 500 and a system size comparable to experimental setup behaves similar to experiments and simulations.Fit of Eq. ( 1) with τ N = 26:2 ± 0:2 min.Details in SI.Source data are provided as a Source Data file.
(initially) dominant.We then find that the time evolution of the average number of dimerized particles can be described by a Poisson process 〈N c (t)〉 = 2rt with a time-dependent rate r = r(t).Assuming a linear decrease of non-dimerized particles, N f = N − 2rt, we find the Michaelis-Menten equation where τ N is the time required for two particles to collide and interlock (see SI).At short times 〈N c 〉/N ~t/τ N increases linearly with t, as the abundance of non-dimerized crescents makes the rate approximately time-independent.At long times, on the other hand, all particles are dimerized and 〈N c 〉/N → 1. Fitting Eq. (1) to our experimental and numerical data we find very good agreement, as confirmed by Fig. 2f, g and Fig. 3b.In the following we will return to using the model of ref. 37  for simulations.The simulations presented are performed in the limit where translational and rotational noise are vanishing because when the rotational noise is small, the change of number of clustered particles is also relatively small.Only for large values the number decreases significantly, see SI for the effect of noise on the clustering dynamics.

Density dependence of crescent clustering
Our initial results demonstrated that significant clustering already appears at very low particle densities.To investigate the density dependence of clustering and explore whether there is a limit, we tested three even lower particle concentrations in experiments.
Remarkably, even at concentrations as low as nine crescents/mm 2 (i.e.ϕ s = 0.022%) cluster formation still happens, see Fig. 3a, albeit slower.This surface area fraction is two orders of magnitude smaller than for active spheres 1,2,17 , rods 29 , and standing disks 32 .We return to simulations to investigate even lower particle densities.Similar to the experiments, a significant number of particles cluster (Fig. S11).Plotting N c /N in Fig. 3b, d, we find that the assembly occurs slower at lower particle densities, but still looks qualitatively similar.The Michaelis-Menten Eq. ( 1) fits all the simulation and experimental data well, independent of the density.Eq. ( 1) suggests the presence of only a single time scale, the constant τ N , for cluster formation.This is surprising since the active Brownian motion of isolated particles, as well as the interaction between particles could, in principle, give rise to other time scales.To explain this we assume that particle collision occurs predominantly in the diffusive regime, where the particles' mean squared displacement is given by hΔrðtÞi 2 = 4Dt, with D the diffusion coefficient.Thus, the  1).The observation time needed to cover the entire sample reflects as an x-error for each data point (also applicable to (c)).c Master curve obtained from (b) by rescaling time as t → tN(t).Fit with Eq. ( 1) with τ 1 = ð2:9 ± 0:3Þ10 4 min.d N c /N(t) from simulations for various densities, and (e) the corresponding master curve, where τ 1 = ð2:0 ± 0:1Þ10 4 min.Source data are provided as a Source Data file.
time for two particles to collide is approximately τ N ∼ ' 2 =ð4DÞ, with ℓ the average inter-particle distance.If particles are uniformly distributed in space, (D/ℓ) 2 ≈ N, with D the system size, thus τ N ≈D 2 =ð4DNÞ.Notice that, consistently with our combinatoric calculation, τ N ~1/N (see SI). Transforming t → tN is then equivalent to rescaling time by the only time scale of the process, thus removing any particle-density dependence.
Indeed, with this simple rescaling of time, we obtain a master curve for the data from experiments and simulations described by a single fit parameter τ 1 , shown in Fig. 3c, e.We attribute the difference in fit parameter τ 1 = Nτ N of approximately an order of magnitude to flowinduced attractions present in experiments as well as to slow sedimentation of particles which increases the surface concentration in time, see SI for a discussion.Both effects only modify the time scale but are not relevant for the basic clustering dynamics.We can thus conclude that the dimerization of active crescents is a single-time-scale process, regardless of the particle density.Furthermore, clustering occurs even at arbitrarily low concentrations.This is in stark contrast to spheres, which undergo a dynamic clustering process that consists of break up and growth, and the transient pairs formed by rods and disks.This peculiarity of active crescents ultimately determines the high performance of their assembly and originates from the efficiency of the interlocking mechanism governing dimerization.This implies an infinite interaction time and, effectively, a vanishing ballistic-diffusive crossover time thus leaving τ N as the only finite time scale of the process.

Shape dependence of crescent clustering
To study the connection between shape and collective behavior we print two other types of bent rods, one with a smaller opening angle α = 90°and one with a larger opening angle α = 260°, keeping L constant, see Fig. 4c and e.We measure their clustering behavior for different surface area fractions and, strikingly, find in both cases that the N c /N data can again be rescaled into a single master curve.The fit parameter τ 1 increases from τ 1 = ð2:9 ± 0:3Þ10 4 min for α = 180°( Fig 3c) to τ 1 = ð3:6 ± 0:3Þ10 4 min for α = 90°(Fig.4d) and τ 1 = ð7:0 ± 0:5Þ10 4 min for α = 260°(Fig.4f), reflecting the slower assembly process.
To better compare these numbers, we also need to take into account the particles' speed, as a higher speed will lead to faster assembly.At 1% H 2 O 2 the speeds for 90°crescents at 〈v〉 90 = 1.02 ± 0.03 μm s −1 were comparable to that of the 260°crescents, 〈v〉 260 = 1.09 ± 0.01 μm s −1 , implying that the 260°crescents assemble significantly slower.However, in both the 90°crescents and the 260°crescent case, the speed is higher than the speed of the 180°p articles, 〈v〉 180 = 0.78 ± 0.08 μm s −1 making the difference in the characteristic assembly time τ 1 even more noteworthy.Moreover, the number density of particles in the sample is higher at the same ϕ s for the 90°crescents, due to their shorter arc-length, underlining the difference further.
To investigate if interlocking due to the concave shape is important, we execute two control experiments.First, we print straight rods, i.e. α = 0°, and quantify their clustering behavior in time for different surface area fractions, see Fig. 4a, b.Straight rods can form pairs too 29 but can still slide past each other leading to a breakup of the cluster Fig. 4i and Supporting Video 9.At 1% H 2 O 2 , the rods move with an average speed of 〈v〉 0 = 1.58 ± 0.03 μm s −1 and, surprisingly, follow the same clustering dynamics as the bent rods with a fit parameter for the master curve τ 1 = ð4:1 ± 0:2Þ10 4 min.This means that the clustering behaviour of active straight rods is within error comparable to the one of 90°bent rods, noting that the speed of the straight rods was higher than of any of the crescents.
The second control group are spherical particles, which are at the other end of the shape spectrum.This shape change can be effectively implemented by crescent particles that move in the direction of and interact with their convex side instead of their concave side, see Fig. 4g, h, with effectively a collision interaction comparable to those of spheres in 2D experiments.We purposefully did not print bent rods upside down to obtain convex-side leading swimmers because such particles would show a truncated shape with a flat surface at their convex side due to the contact area between printed structure and substrate required in the printing process.This flat surface leads to a significantly higher stability of pairs and thus affects the clustering dynamics.Instead, we experimentally achieve convex-side leading crescents by increasing the H 2 O 2 concentration to 5%, which leads to a direction reversal, and term these c + crescents, following ref.27.While their average velocity at 1.3 ± 0.2 μm s −1 was slightly higher than that of the concaveside-leading particles, their efficiency to cluster at low density drops significantly (Fig. 4h).The few pairs observed over the course of the experiment typically broke up by sliding past each other within ~2 min (Fig. 4i and Supplementary Video 12).Fig. 4i shows pairs of straight and bend rods of different opening angles and which all differ in their ability to interlock (also see Supplementary Video 7-11).For all particle shapes that significantly slow down the motion, i.e. straight rods and concave-side leading bent rods, we observe clustering despite the low surface densities.This motility induced clustering results from a combination of the break down of detailed balance and steric interactions.In bulk systems, a similar mechanism drives phase-separation at small length scales (MIPS).
Here we showed that the same mechanism can lead to a cluster phase even in the absence of a true bulk phase, and hence term it motility induced clustering.
To study many different opening angles we turn again to simulations.We consider eleven different opening angles from 30 °to 330°i n steps of Δα = 30°, keeping the cross-section fixed, see Fig. 5a.For fixed arc-length and varying cross-section, the particles cluster at a similar speed, see SI.When again plotting N c /N we find that 180°c rescents cluster most efficiently (cf.Fig. 5b) with the smallest value for τ 1 .To quantify this behavior we present in Fig. 5c the inverse fit parameter 1/τ 1 as a function of α, normalized to the value found for α = 180°.The larger 1/τ 1 , the faster the clustering.We find the largest value for 180°crescents, with 1/τ 1 being approximately symmetric about α = 180°up to 90°.However, increasing the opening angle further, τ 1 increases sharply, while with decreasing opening angle the fit parameter reaches a plateau around α = 90°.
To rationalize the dependence of τ 1 on α we assume two competing effects are relevant.This approach is similar to ref. 43 where, however, a three-dimensional system of passive bent rods was considered which results in significantly different dynamics.1) the probability P cl to form a cluster and 2) the stability of the cluster.P cl depends on the opening length L op , which is independent of α for α < 180°.However, for α > 180°, L op and thus P cl , decrease with increasing α.In the extreme example of a sphere, stable dimers cannot form.The probability of two particles colliding at an angle such that their cavities are facing each other, resulting in a dimer, is However, once a cluster has formed, its stability depends on α, too.For smaller α, particles interlock less and cluster are less stable, q ÞΘðπ À αÞ + 0:106P cl Θðα À πÞ, with Θ the Heaviside step function.The red data points (star) represents the normalised inverse fit parameter for straight rods, 90 ∘ crescents and 260 ∘ crescents (from left to right) found in experiments before rescaling, the blue data points (diamond) the value after rescaling (see main text and SI).c Evolution over time of the number of particles in a two-particle-cluster (top), a three-particle-cluster (center) and a fourparticle-cluster (bottom), respectively.Source data are provided as a Source Data file.
e.g., in case of a collision with a free particle.As a simple proxy for stability we define the length dðαÞ = ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi ðL op =2Þ 2 + h 2 q that one crescent in a dimer would have to move to leave the cluster.The smaller this length, the more unstable the cluster.This length scale yields a time scale T ∼ d=v 0 such that τ 1 ∼ T captures the α-dependence of the parameter τ 1 .We find 1/τ 1 ~1/d(α).We find that combining these two effects results in a good agreement, see Fig. 5c, especially considering the simplified arguments we used.
Cluster formation is optimal for 180°crescents because the probability for cluster formation is high, due to their relatively large opening length, and, once formed, clusters exhibit significant stability from the particles' strong curvature.In experiments we also find 1/τ 1 to be highest for 180°particles.However, the values found for straight rods, 90°crescents and 260°crescents are significantly higher than the values found in simulations, see blue data points in Fig. 5b.We can rationalize this difference from the different propulsion speeds and sedimentation rates of those particles.By rescaling the inverse and normalized fit parameters τ 180 1 =τ 1 with the measured propulsion velocities for the respective shapes as well as with the corresponding inverse calculated sedimentation velocities (see SI and 44 ), we obtain corrected values for τ 180 1 =τ 1 which are in perfect agreement with the values found in simulations, implying again that hydrodynamic interactions are negligible, see red diamonds in Fig. 5b.
As a consequence of the shape-dependent competition between interlocking and breaking-up, higher-order cluster are also more stable for 180°.Three-particle-clusters are already most prevalent for 180°p articles, but four-particle-clusters are effectively only observed for α = 180°(see Fig. 5c).While the absence of higher order clusters in samples with higher opening angles can be explained by the increased difficulty to form these types of structures, for lower opening angles, the lack of larger clusters arises from the instability of two-and threeparticle clusters when colliding with free particles.

Discussion
In this manuscript we investigated the clustering behavior of selfpropelled particles with a shape that is continuously tunable along a single anisotropy dimension by exploiting a combination of 3D printing and simulations.Printing bent rods with different opening angles allowed connecting our results with the two major shapes thus far employed, spheres and rods.Surprisingly, stable interlocking into pairs occurs at arbitrarily low number concentrations of particles due to a interlocking type interaction.Pair formation plays a critical role in nucleation, and the dynamics of the fraction of clustered particles can thus be captured with a notably simple formula which allows rescaling of the clustering dynamics into a single master curve.We find that particle shape influences both the assembly speed and cluster size distribution with an optimum for crescents with an 180°opening angle and can be described by a competition between interlocking and break-up.
Our experimental and numerical results provide the first detailed understanding on how a gradual shape change from a sphere to a rod affects the clustering behavior of active particles and uncover a simple yet general description for this complex effect.These key insights are important to understand how non-spherical shapes affect biological swimmers, and can be leveraged to further design interlocking type connections between active particles to control the stability of complex self-assembling structures through the depth of their entangling sites independent of their composition and surface chemistry.This strategy offers the unprecedented opportunity to design hierarchical assembly pathways towards active functional materials through engineering the active particle shape.The guaranteed assembly even at extremely low concentrations together with the interlocking shape recognition furthermore opens up possibilities in the drug delivery where site recognition or a low concentration of the drug-carriers as well as the activation of reactants after assembly is important.In the future, it will be interesting to study the cross-over from our lowdensity system which exhibits motility induced clustering of mostly small clusters to a system with higher surface area density where larger clusters with different assembly dynamics and motility induced phase separation (MIPS) may occur.

Reagents
Fused silica substrates as well as the photoresist IP-Dip were acquired from Nanoscribe GmbH.Propylene glycol methylether acrylate (PGMEA, ReagentPlus ≥99.5%) was purchased from Sigma Aldrich and isopropanol (IPA) was obtained from VWR chemicals.Hydrogen peroxide (H 2 O 2 , 35 wt% solution in water, stabilized) was purchased from Acros Organics.Water was purified by means of a MilliQ system (resistivity ≥ 18 MΩ.cm).Unless otherwise noted, all chemicals were of analytical or reagent grade purity and were used as received from commercial sources.

Particle fabrication
All particles were 3D printed on a Nanoscribe Photonic Professional GT microprinter equipped with a 63x oil-immersion objective (Zeiss, NA = 1.4) following 35 .Particles were printed on a fused-silica substrate using IP-Dip as photoresist.Bent rods with an opening angle of 180°a nd 260°were designed with a thickness of 2 μm and a cross-section length L of 10 μm.For bent rods with an opening angle of 90°the thickness of the particles was reduced to 1.5 μm in order to keep the same size ratio between the thickness of the crescent legs and the concave opening.Straight rods were designed with a thickness of 2 μm and a length L of 10 μm.During development the printed particles were submerged in PGMEA for 30 min followed by 2 min in IPA and then left to dry overnight at ambient conditions.Once dry, active crescents were coated with a 5 nm layer of catalyst (Pt/Pd 80:20) using a Cressington 208HR sputter-coater.To prevent excess Pt/Pd, the area around the print was protected with tape which could be removed later without damaging the printed structure.After sputter coating, each print was checked by default under the microscope to determine its quality.The substrate was placed onto a 200 μL MilliQ water droplet in the center of a small glass petri dish such that the printed crescents were in contact with the water.Repeated sonication for 2 min allowed recovery of 90% of the particles from the substrate.Particles were subsequently concentrated by centrifugation followed by removal of the supernatant.

Imaging, data acquisition, and data analysis
For observation of the clustering dynamics, a Nikon Eclipse Ti-E bright-field light microscope with a Plan Apo λ 20x long working distance objective (NA = 0.75) was used.Pt-coated colloidal particles were suspended in freshly prepared 1% hydrogen peroxide aqueous solution in the case of concave-side-leading crescents and rods, or in 5% hydrogen peroxide for convex-side-leading crescents.Control experiments with Brownian crescents were performed in MilliQ water.This colloidal solution was then placed into a sample holder (+ = 8 mm) using an untreated borosilicate glass coverslip (VWR, 25 mm, No. 1) as substrate.
The cluster distribution in the sample was recorded after different time steps.We hereby defined a cluster as an assembly of at least two particles that remained stable for at least 5s.A square array of 7 × 7 images (2048 × 2048 px) with an overlap of 1% was taken to cover and observation area of 4.5 mm x 4.5 mm in the central area of the sample holder.The time needed to take these 49 single images including corresponding 5s videos for each field of view resulted in an average value and error for the determination of the time-point.Short videos were only taken for samples with lower particle densities where the total observation time exceeded 100 min.All measurements were taken in the dark and in the absence of bubble formation.
Single crescents as well as crescent cluster species were counted for every time step, where single stuck particles (as determined from 5s videos) were left out of the count of single crescents.Average particle velocities were captured from 30s long videos with a frame rate of 20 fps.Particles were tracked in each frame, applying the Canny Edge detection algorithm to generate a mask out of which the particle center of mass was obtained.Individual crescent-velocities v were determined by fitting the short-time regime (Δt ≪ τ r ) of the mean squared displacement (MSD) of different crescents with Δr 2 = 4DΔt + v 2 Δt 2 .The second fit parameter is the diffusion coefficient D and τ r was chosen as the rotational diffusion time for a sphere.We here use the aforementioned fit-equation, originally developed for spherical active particles 36,45 as an approximation for our crescentshaped particles despite them not being spherical.

Dynamic simulations
Following Wensink et al. 37 we simulate N self-propelling particles in two dimensions that move with a velocity v 0 .Each particle is discretized into k equidistant spherical segments of diameter d (see Supplemental Figure 2).A particle ρ has an orientation u ρ and position r ρ .The position of segment i with respect to the position of the center of mass r ρ is denoted by e i ρ .The pair potential of two particles ρ and δ is given by U ρδ = k À2 P k i,j = 1 uðr ij ρδ =dÞ.Here, uðxÞ = u 0 exp Àx ð Þ=x 2 is a short-range potential that is repulsive if u 0 > 0 resulting in effectively hard particles.r ij ρδ = |r ρ À r δ + e i ρ À e j δ | is the distance of two segments of the two different particles.The equations of motion, found from balancing forces and torques due to activity and steric interactions, for a particle ρ are then given by 37 : where f t and f r are translational and rotational friction, respectively, u ρ = fsin φ ρ , cos φ ρ g, v 0 = F a /f t , and U = ∑ ρ,δ(ρ≠δ) U ρδ /2.ξ and η are translational and rotational Brownian noise, respectively.We simulate the particles in a square system with periodic boundary conditions.A number of N particles is initialized in a system of size D × D with random orientation and position.In the simulations we have defined a particle to be in a cluster with another particle if it is closer than 3R for 100 iterations.Unless otherwise noted, each of the results presented is found by averaging over 100 independent runs.We employed the same particle dimensions and velocity in the simulations as in experiments.A list of parameters used, the mapping between computational and physical units, as well as additional details can be found in the Supplementary Methods section.

Balls-into-bins model
The dynamics of N active particles of cross-sectional length L moving with speed v 0 in a square system of width D are mapped onto a ballsinto-bins model.For this, we choose the number of boxes to be B = (D/S) 2 , with S ≈ L the size of a single box.In the first step of the simulation we randomly assign each of the N balls a box.If n > 1 balls are assigned to the same box, we count this as a n-cluster.If a ball is in a box by itself, it is counted as a free particle.In the following step of the simulation the above procedure is repeated, but only free particles are assigned a new box while n-clusters remain in their respective box.See the Supplementary Methods for a more detailed explanation, a list of parameters used, and a mapping between computational and physical units.

Fig. 1 |
Fig. 1 | Experimental design of self-propelled bent rods.a By changing the design parameters, the shape of the anisotropic microswimmers can be continuously tuned from a sphere to a bent or straight rod.Cross-sectional length L, arc-length L arc , opening angle α, and radius of curvature R. b Schematic illustration of 3D printing of a bent rod by two-photon polymerization.c Scanning electron microscopy image (top left inset) and bright-field image of a 3D printed active crescent with α = 180 ∘ , see the Methods section for details on particle imaging.The particle moves, concave-side leading, in the direction of the red arrow.Red line is the trajectory of past 30 s.

Fig. 2 |
Fig. 2 | Time-dependent clustering of active crescents.a Bright-field microscopy image taken 90 min after mixing particles and fuel solution shows that interlocked pairs and clusters, as well as single particles are present (surface-area fraction ϕ s = 0.261%, corresponding to 104 crescents/mm 2 ).b Simulation snapshots and bright-field microscopy images of clusters of different size N. c Snapshot of simulations for ϕ s = 0.261%.Clusters are colored according to size.d Image sequence of formation of a 180 ∘ crescent-pair in time.The pair only rotates after interlocking.e Time evolution of the probability P(n) to find a crescent in a cluster of size n

Fig. 3 |
Fig. 3 | Density dependent clustering of concave-side leading, bent rods (α = 180 ∘ , L = 10 μm). a Cluster-size distribution over time for different surface-area fractions ϕ s , corresponding to particle number densities of (from left to right) 9, 15, and 25 crescents/mm 2 .Note that ϕ s is time-dependent and the stated values are measured after 90 min (see main text and SI).(b) N c /N(t) for different densities fitted with Eq. (1).The observation time needed to cover the entire sample reflects as an x-error for each data point (also applicable to (c)).c Master curve obtained from (b) by rescaling time as t → tN(t).Fit with Eq. (1) with τ 1 = ð2:9 ± 0:3Þ10 4 min.d N c /N(t) from simulations for various densities, and (e) the corresponding master curve, where τ 1 = ð2:0 ± 0:1Þ10 4 min.Source data are provided as a Source Data file.

Fig. 4 |
Fig. 4 | Shape-dependent clustering behavior.Scanning electron microscopy (inset) and bright-field images of 3D printed (a) straight rods (α = 0°) and crescents with (c) α = 90°, (e) α = 260°and (g) α = 180°.The trajectory travelled by the active particle in the past 30 s is shown in red and the direction of motion is indicated by a red arrow.Fraction of particles in a cluster plotted after rescaling time as t → tN(t) for systems of (b) self-propelled straight rods, (d, f) concave-side leading crescents with (d) α = 90°and (f) α = 260°, and (h) convex-side leading crescents with α = 180°.The fit parameter obtained from fitting with Eq. (1) is (b) τ 1 = (4.1 ± 0.2)10 4 min, (d) τ 1 = (3.6 ± 0.3)10 4 min, and (f) τ 1 = (7.0 ± 0.5)10 4 min respectively.The xerror indicates the observation time needed to measure the entire sample.All surface area fractions ϕ s correspond to particle densities below 50 crescents/mm 2 .The exact values can be found in the SI.i Bright-field images of rods and crescent pairs that differ in their ability to interlock, cf.Supporting Video 7-12.The white arrow indicates how pairs that are not sufficiently stabilized by interlocking can slide past each other and break.Source data are provided as a Source Data file.

Fig. 5 |
Fig. 5 | Simulation results on the influence of the opening angle α on the clustering behaviour.a Bent rods with increasing opening angle from 30 ∘ to 330 ∘ used in the simulations shown in (b).b Inverse fit parameter 1/τ 1 as a function of the opening angle α, normalized to the value of τ 1 for α = 180 ∘ .The fit is found from a simple model describing cluster creation and decay, see the main text for details.The fit shown here is τ 180 1 =τ 1 = 2:357ð1:149 À 1= ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi dðαÞ 2 + 1